Universal corrections to the Fermi-liquid theory 
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We show that the singularities in the dynamical bosonic response functions of a generic 2D Fermi 
liquid give rise to universal, non-analytic corrections to the Fermi-liquid theory. These corrections 
yield a term in the specific heat, T terms in the effective mass and the uniform spin susceptibility 
Xs{Q = 0,T), and \Q\ term in Xs(Q,T = 0). The existence of these terms has been the subject of 
recent controversy, which is resolved in this paper. We present exact expressions for all non-analytic 
terms to second order in a generic interaction U{q) and show that only (7(0) and U(2pf) matter. 

PACS numbers: Tl.lOAy, 71.10Pm 



The universal features of a Fermi liquid and their phys- 
ical consequences continue to attract the attention of the 
condensed-matter community. In a generic Fermi liq- 
uid, the imaginary part of the retarded fermionic self- 
energy Y,ii{k,uj) on the mass shell is determined solely 
by fermions in a narrow uj) energy range around the 
Fermi surface and behaves as S" oc -|- (ttT)^ P|. This 
regular behavior of the self-energy has a profound effect 
on such observables as the specific heat and uniform spin- 
and charge susceptibilities, which have the same func- 
tional dependences as for free fermions, i.e., the specific 
heat is linear in T and the susceptibilities approach fi- 
nite values at r = 0. A regular behavior of the fermionic 
self-energy is also in line with a general reasoning that 
turning on the interaction in I? > 1 should not affect 
drastically the low-energy properties of a system 3| , un- 
less special circumstances, e.g., a proximity to a quantum 
phase transition interfere. 

The subject of this paper is the analysis of non-analytic 
corrections to the Fermi- liquid behavior in a generic, 
clean Fermi liquid. These corrections are universal in 
a sense that they are determined by fermions near the 
Fermi surface. It has been known for some time that cor- 
rections to the Fermi-liquid form of S"(w) do not form a 
regular, analytic series in but rather scale as ui^ for 
2 < D < 3, with an extra logarithm in Z) = 2 [1 H |1| . 
[For \ < D <2 this form persists, but it is not a "correc- 
tion" anymore.] These non-analytic u)^ terms (as well as 
non-analytic vertex corrections) are of fundamental in- 
terest as they may give rise to anomalous temperature 
and momentum dependences of observable quantities. A 
well-known example is the In T term in the specific 
heat of a 3D Fermi liquid, caused by the anomalous term 
S"(u;) oc oj^ (and hence S'(a;) ex w^logw) A related 
example is the linear-in-T correction to the conductivity 
of a weakly disordered 2D system 0, S 01 ■ Non-analytic 
corrections are also important for the theory of quantum 
critical phenomena in itinerant ferromagnets j^], as a 
non-analyticity of the static spin susceptibility changes 
the nature of the phase transition . 

Belitz, Kirkpatrick and Vojta (BKV) fT5| and, later. 



Misawa |0| argued that the non-analyticity in the 
fermionic self-energy should gives rise to a non-analytic 
momentum expansion of the particle-hole susceptibility 
x{Q,T). For non-interacting fermions, x(Q,0) is given 
by the Lindhard function and is analytic in Q for small 
Q in all D. Diagrammatically, corrections to the Lind- 
hard function are obtained by self-energy and vertex- 
corrections insertions into the particle-hole bubble (cf. 
FigI2l. Diagrams with self-energy insertions are readily 
estimated by power counting, and the result is that the 
non-analytic, uj^ term in E gives rise to 5x{Q) Q^~^ 
,with extra logs for D = 3 and D = 1. Power count- 
ing also suggests |0, that the non-analyticity in S 
should affect the temperature dependence of the uniform 
susceptibility x(0,r) and gives rise to Sx{0,T) cx T^^^, 
with extra logarithms in D = 3 and D = 1. By the same 
arguments, non-analyticity in E should lead to the T^~^- 
dependence of the effective mass m*{T) and to the T^- 
dependence of the subleading term in the specific heat, 
SC{T), with extra logarithms in 3D and ID 14, 15]. 

Our motivation to pursue a further study of non- 
analytic corrections to the Fermi liquid behavior is two- 
fold. First, we want to verify that power counting argu- 
ments by carrying out an explicit analytic calculations 
of several observable quantities: specific heat, effective 
mass, and spin- and charge susceptibilities. That power 
counting may be misleading is seen, e.g., from the exam- 
ple of the free-fermion susceptibility: according to power 
counting, it should also have a non-analytic momentum 
dependence, whereas the exact result is analytic in Q for 
small Q. The existing literature on this issue is contro- 
versial. BKV verified their power counting arguments 
for the spin susceptibility in 3D by explicitly comput- 
ing (5xs(Q,T) to second order in the interaction. They 
demonstrated that Sxs{Q,0) oc Q^lnjQ], in agreement 
with power counting. At the same time, Carneiro and 
Pethick 0] and later BKV found that uniform Xs(0,T) 
scales as but not as T^lnT, as predicted by power 
counting. On the contrary, Misawa did find a T^lnT 
term in his calculation [Tg]. In 2D, BKV conjectured 
that Xs((5,0) scales as ]Q] but no explicit calculation has 



not been performed. Hirashima and Takahashi com- 
puted x(0, T) in 2D numerically, but coud not draw any 
definite conclusions about the T-dependence because of 
numer ical difficulties. Chitov and MiUis (CM) [13, found 
analytically that the leading term in Xs (0, T) in 2D scales 
as T, in agreement with power counting [l3j. At the same 
time, CM found that different contributions to the non- 
analytic T term in m*{T) and to the term in 5C{T) 
(both predicted by power counting), cancel each other, 
and only analytic corrections survive. Meanwhile, Bedell 
and Coffey [l5j reported a term in SC{T). Very re- 
cently. Das Sarma et al. found a linear-in-T term in the 
effective mass for the Coulomb interaction in D = 2 . 

In this paper, we present analytic results for the spe- 
cific heat, effective mass, and spin and charge suscep- 
tibilities of an interacting 2D Fermi system, up to sec- 
ond order in the short-range interaction U{q). We found 
that power counting arguments are generally valid, i.e., 
Xs{Q) « IQI, XsiT) oc T, m(r) oc 7\ and 5C(r) cc T\ 
These results agree with Refs. [IJ lli [IJ; the form of 
Xs{T) agrees with that found by CM but the forms of 
miT) and 5C{T) disagree with those by CM. Still, our 
prefactors for non-analytic terms differ from that found 
in [H [13. We also verified that in 3D, Xs(Q,0) oc 
Q^lnQ while Xs(0, T) oc T^, in agreement with BKV. 
Finally, in agreement with CM, we found no non-analytic 
terms in the charge susceptibility, Xc- 

Another motivation for this study is to clarify the ori- 
gin of the non-analytic corrections in the Fermi-liquid 
theory. We found that these corrections originate from 
the singularities in the dynamic particle-hole response 
function, Tl{q, 51), near q = Q and q = 2fci?, where n(g, 51) 
is non-analytic. For D = 2, near q = 
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where q = q — 2kF and q ^ 2kF. Physically, these two 
singularities give rise to a zero-sound mode and Friedel 
oscillations, respectively. The singularity near g = is 
entirely dynamic, while the one near 2kF is also present 
in the static limit for q > 2kF- We found that the singu- 
larities in n((7, 51) are necessary ingredients which make 
power counting arguments valid. Furthermore, we found 
that the singular pieces in the effective mass, specific 
heat and Xs{Q,T) originate exclusively from the scat- 
tering amplitude with zero momentum transfer and zero 
total momentum: Ta.p-'y.sik, —k; k, —k) = U{0)5af5j3s — 




b) 




k p p+q k+q k 



FIG. 1: Non-trivial second-order diagrams for the self-energy. 



U{2kF)SasSf3'y This implies that (i) non-analytic terms 
depend only on U{Q) and U{2kF) but not on the interac- 
tion averaged over the Fermi surface, and (ii) up to over- 
all sign, the non-analyticities at g « and q k 2kF con- 
tribute equally to individual diagrams for the fermionic 
self-energy and Xs, i-e., singular corrections to a Fermi 
liquid can be viewed equivalently as coming either from 
the singularity in the dynamical particle-hole bubble at 
q — or at q = 2kF ■ As the q = singularity is entirely 
dynamical, the non-analytic corrections to a Fermi liquid 
are dynamical in nature as well. 

Effective mass and specific heat. To find the ef- 
fective mass, m*{T), and the correction to the specific 
heat, dC{T), one needs to know the real part of the 
fermionic self-energy, S'(fc,Ci;), on the mass shell, i.e., at 
(k^ — kp) /2m = Lo. The two nontrivial second- 



er 



order diagrams for the Fermi energy are presented in 
Fig. n We evaluated E"(fc,Ci;) first, using the spec- 
tral representation, and then obtained S'(a;) on the mass 
shell via Kramers-Kronig transformation. The imaginary 
part of the self-energy reduces to well-known forms [5j 
Y."(k,uj) oc intj and S"(fc,cj) oc InT for k near the 
Fermi surface and in the limits of T — > and w ^ 0, 
respectively. We obtained Yj"{k,Ljj) at arbitrary w/T and 
Lo/tk- The full expressions are, however, rather involved 
and will be presented elsewhere [2^. Using full S"(fc, w), 
we find for the real part of the self-energy on the mass 
shell, neglecting a regular, Fermi-liquid type lo term 



S'(c.) 



mU 



IGtt' 



where = U^{Q) + U^i2pF) - U{0)U{2pf) and 



g{x) = 1 + 4 



^+Li.(-e-l^ 



(3) 



(4) 



where Li2(a;) is a polylogarithmic function. 

In the two limits, g{oo) = 1 and g{x ^ 1) « 41n2/a;. 
The first limit corresponds to T = in which case Eq. 
gives S'(w) oc uj\uj\. This non-analytic form agrees 
with power counting. For small lu/T, i.e., for a: <C 1, 
the 1/x form of g{x) leads to the ojT term in S'(w) for 
LO <^T . This in turn implies that the quasiparticle mass 
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m*(T) acquires a linear-in-T correction 



In 2 T 



(5) 
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Using Eas. pHl . we find a correction to the specific heat 



(5C(T) 



2m 



ouj 



= 0.f74 Cfl m^U^ 



EfJ 
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where Cfl = ttTto/S is the Fermi-gas result. We see 
that a non-analyticity in the fermionic self-energy gives 
rise to the term in the specific heat. This term comes 
only from fermions in a near vicinity of the Fermi surface 
and from this perspective is model-independent. 

We also verified that the non-analytic part of 
stems exclusively from the scattering process in which 
two internal momenta in the self-energy diagram are near 
— k, whereas the third one is near k, i.e. , when both 
the total and transferred momentum are near zero. As 
an independent check, we obtained the non-analytic part 
of ^'{(jj) by re-expressing the self-energy in terms of the 
particle-particle bubble-Eq. then results from the 
(logarithmic) singularity of the particle-particle bubble 
at small total momentum and frequency. This should 
indeed be the case if only the scattering amplitude with 
zero total momentum is important. 

Spin susceptibility. The relevant diagrams for the spin 
susceptibility are presented in Fig|21 Evaluation of the 
diagrams is rather tedious but straightforward. We cal- 
culated all diagrams in two ways: i) explicitly, by ex- 
ploring the non-analyticities in the particle-hole bubble 
near q = and 2kF, and ii) by retaining only vertices in 
which both total and transferred are small. In the sec- 
ond approach, we expanded in total and transferred mo- 
menta and extracted universal contributions, which are 
independent of the upper cutoff for the expansion. We 
obtained identical results in both methods, which proves 
that only a single scattering amplitude is relevant. 

The non-analytic contributions to the spin susceptibil- 
ity from individual diagrams are as follows 



XiiQ,T) = xoK{Q,T)[U^iO) + U^2kF)] 

X2iQ,T) = -xoK{Q,T)U{0)U{2kFy, 

X3{Q,T) = xoKiQ,T)[U^i2kF)-U^iO)] 

X4{Q,T) = xoK{Q,T)UiO)Ui2kFy, 

X5 = X6 = X7 = 0, 



(7) 



where xo = ttl/tt, and K{Q, 0) and K{0, T) are given by 
2 fmUV \Q\ , , fmUV T 

(8) 










FIG. 2: Relevant second-order diagrams for spin- and charge 
susceptibilities. The last two diagrams are non-zeroes only 
for the charge susceptibility. 



Collecting all contributions, we find 
7 

Xs{Q, T) = Y, x^{Q, T) = xo [1 + 2K{Q, T)U\2kF)\ . 

i=l 

(9) 

We see that all non-analytic contributions with U{0) can- 
cel out, and the final result depends only on [/(2A;_f). 

We also performed a similar calculation in 3D and re- 
produced the BKV's result-the analog of Eq. jnj but with 
K3DiQ,0) = (l/18)(mfcFt//47r2)2(Q/fc^)2ln(A:j./Q). In 
agreement with BKV, we also found that K{0,T) oc 
with no logarithmic corrections. 

We now look more deeply into how the non-analytic 
contributions to susceptibility emerge. The power- 
counting argument does not rely on the singularity in 
the particle-hole bubble. Indeed, near q = 0, the singular 
f^rn/V {vFqy + piecc in Tiph[q,^m) has the scaling 
dimension of one and hence can be treated as a constant 
in power counting. However, we found that for each di- 
agram, a replacement of T\ph{q,^lm) by a constant does 
not give rise to a linear-in-|(5| term in Xs('5,0) because 
the prefactor for the \Q\ term contains the integral over 
g, in which all poles are located in the same half-plane. 
The q integral then obviously vanishes. This vanishing 
could not be detected in power counting. The substitu- 



4 



tion of the full Ilph{q,i^m) into the susceptibility makes 
power counting arguments valid as Ilph{q-,i^m) contains 
a branch-cut singularity which is present in both half- 
planes of q. With this term present, the location of the 
poles in a complex q plane becomes unessential, and the 
q integral does not vanish. 

The linear-in-T dependence of x(0, T) is also associ- 
ated with the singularity in Ilph{q,^), but the way it 
emerges is different from T ~ and is similar to an 
anomaly. Consider for example the q — contribution to 
Xi(0,T). Leaving the integration over q and summation 
over flm as the last operations, we obtain 



=°(0,T) = -2xo 



mt/(0) y 
47r J 



where 



dqvpq 



nl{2nl-{vFqr) 



{vpq) 



2« 



(10) 



(11) 



Contrary to the T = case, the result for J at finite 
T depends on the order of the integration and summa- 
tion. Indeed, integrating first over q in infinite limits one 
obtains 1/4 for all values of and a subsequent fre- 
quency summation does not yield a universal piece con- 
fined to low energies. On the other hand, performing the 
summation over f2m first, keeping q finite, and then in- 
tegrating over g, one obtains a universal —1/4 piece in J 
which gives rise to a linear-in-T piece in Xs(0,T). The 
correct way to compute J is the second one, because the 
f2m-independent result of the integration over q is incon- 
sistent with the fact that the integrand in (|ll|l obviously 
vanishes at f2„i = for all finite q. In reality, q is always 
bounded from below by the inverse system size. 

We see that despite a formal analogy between the 
Q- and T- dependences of Xs^ the mechanism behind 
Xs(0,T) oc r is very different from the one that leads 
to Xs{Qi 0) c)c IQI as in the latter case the order of inte- 
gration is unessential. The factor J in 1)1 0|l can be readily 
found for arbitrary D. For D > 2, we have 



J 



{D-2){A-D) 
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T 



D-2 
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D-2 



1 



(12) 



For _D ^ 2, we reproduce J — > —1/4. For arbitrary 
2 < £1 < 3, J oc T^-^, i.e., Xs {T) oc T^^^. For D = 3, 
however, the integral is regular, and J (xT, i.e., Xs{T) oc 

without logarithmic corrections. 
Charge susceptibility. For the charge susceptibility, we 
have two additional contributuions given by diagrams 6 
and 7 in Fig. 2. We find 

X6{Q,T) = -xo K{Q,T) [C/2(0) - C/2(2fcf )] ; 
X7{Q,T) = -xo K{Q,T) [uH0) + U^{2kF)] , (13) 

where K{Q, T) is given by (|SJ|. Combining this last result 
with Eq. Q, we find that all non-analytic terms from 



individual diagrams cancel out, i.e., the charge suscepti- 
bility is regular, in agreement with Ref.|l9|. 

To conclude, in this paper we demonstrated that the 
universal singularities in the bosonic response functions 
of a Fermi liquid give rise to universal non-analytic cor- 
rections to the Fermi-liquid forms of the self-energy and 
thermodynamic variables. We obtained explicit results in 
2D for 5C{T) oc T^, Xs{Q.T = 0) cx |Q|, Xs{Q = 0,r) cx 
T. We demonstrated that these non-analytic terms come 
from the processes with both transferred and total mo- 
mentum close to zero. We also demonstrated that ther- 
mal (oc T)) and quantum (oc |Q|) corrections to the spin 
susceptibility are of different origin. This explains why 
in 3D, Xs((5,0) oc Q^lnQ, while Xs(0,T) ocT. 
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